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Previously  reported  experiments  with  a  self-propelled  body  submerged  in  a  fluid 
with  a  stable  vertical  density  gradient  have  demonstrated  that  the  turbulently 
mixed  wake  first  expands  more  or  less  uniformly  and  then  collapses  vertically 
while  continuing  to  expand  horizontally  (Schooloy  &•  Stewart  1963).  It  was  also 
shown  that  the  vertical  collapse  of  the  wake  generates  internal  waves.  Essen¬ 
tially  two-dimensional  experiments  have  also  been  used  to  explore  some  of  the 
build-up  and  decay  characteristics  of  vertical  wake  collapso  induced  by  a  sub¬ 
merged  burst  of  turbulent  mixing  (Wu  1969;  Schooloy  1968):  The  present  paper 
reports  now  experimental  measurements  and  a  linear  theoretical  analysis  of  tlio 
internal  wave  field  created  in  stratified  water  by  a  burst  of  submerged  turbulent 
mixing.  The  forcing  function  has  beer,  obtained  in  integral  form  for  an  initial- 
value  model  of  wako  collapse  in  terms  of  a  general  Brunt-Vaisiilii  frequency 
profile,  using  normal  mode  theory.  Numerical  results  have  been  determined  for 
the  specialized  case  of  a  completely  mixed  circular  wako  in  a  constant  Brunt- 
Vaisiilu  profile.  These  results  are  compared  to  the  experimental  measurements. 


1.  Introduction 

The  concept  of  internal  wave  generation  by  collapso  of  a  region  of  density 
stratified  fluid  which  has  been  mixed  has  existed  for  many  years  and  various 
researchers  have  contributed  to  the  experimental,  theoretical  or  to  both  sides  of 
the  problem.  Schooloy  &  Stewart  (1963)  measured  the  initial  phase  of  collapso  in 
the  wako  of  a  self-propelled  body  and  demonstrated  that  wave-like  effects 
appeared  at  the  surface.  Furthermore,  they  convincingly  showed  that  the  general 
characteristics  of  tho  surface  phenomena  can  bo  explained  by  the  existence  of  a 
number  of  modes  of  internal  waves  generated  by  tho  wako  collapso.  Their 
analysis,  however,  was  not  aimed  at  predicting  tho  amplitude  of  tho  internal  wave 
motion  from  features  of  tho  initial  mixed  region.  Stockhausen,  Clark  &  Kennedy 
(1966)  also  made  observations  relating  to  wake  collapso  behind  a  self-propelled 
body.  They  wore  concerned  c'niolly  with  tho  shape  of  the  wake  during  tho  initial 
expansion  and  subsequent  collapso  and  did  not  analyse  their  results  from  the 
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point  of  view  of  internal  wave  generation.  Observations  on  the  details  of  a  mixed 
region  during  collapse  in  a  two-dimei«»ional  laboratory  model  were  also  taken  by 
Wu  (I960).,  The  mixed  region  in  this  case  was  formed  by  stirring  a  contained 
semi-cylindrical  region  at  one  end  wall  of  a  tank  containing  stratified  fluid'.  The 
collapse  was  simulated  by  quickly  and  smoothly  removing  the  containing 
apparatus  and  allowing  the  mixed  fluid  to  flow  into  the  stratified:  fluid.  Very 
detailed  results  were  obtained  showing  the  subsequent  shape  of  the  mixed  region. 
The  internal  wave  patterns  associated  with  the  collapse  were  also  determined  and 
results  giving  the  phase  configuration  were  obtained.  No  amplitude  information 
was  illustrated  ^however,  it  was  observed  that  wave  heights  as  large  as  0-2  of  the 
mixed  region  radius  v/erp  formed. 

A  wide  range  of  theoretical  results  have  been  obtained  concerning  internal 
waves  in  general.  Most  results  have  been  obtained  for  linear  models.  A  very  com¬ 
prehensive  treatment  of  the  problem  (for  general  transient  or  steady  generation)  is 
given  by  Lighthill  (1967)  and  further  detailed  theoretical  work  with  careful 
experimental  verification  is  given  in  Mowbray  &  Rarity  (1967).  The  former  is 
concerned  with  a  general  treatise  of  dispersive  waves  in  homogeneous  media  and 
one  section  deals  specifically  with  internal  waves  generated  behind  a  vertically 
moving  steady  disturbance.  Experimental  verification  of  the  phase  configuration 
is  provided.  The  latter  publication  provides  a  complete  linear  theoretical 
treatment  for  forced  internal  waves  in  a  (predominantly)  homogeneous  medium 
and  deals  at  some  length  with  the  asymptotic  solution  of  the  Cauchy-Poisson 
problem.  Tho  observational  data  provided  confirms  the  theoretical  predictions 
regarding  phase  configurations.  No  comparison  of  amplitudes  was  undertaken. 
Impulsively  generated  internal  waves  in  the  atmosphere  have  been  studied  quite 
intensively  (Pierce  &  Posey  1970)  but  not  from  the  point  of  view  of  generation  by 
wake  collapse.  Very  recently,  a  linear  theoretical  analysis  has  been  performed  on 
wave  generation  by  wake  collapse  (and  other  causes)  by  Miles  (1970).  A  non¬ 
linear  numerical  analysis  of  the  problem  has  been  undertaken  by  Wessel  (I960) 
with  a  comparison  to  the  experimental  results  of  Wu  (1969). 

Tho  purposo  of  tho  present  paper  is  to  outline  further  experimental  results 
concerning  waves  generated  by  mixed  region  collapse  and  to  present  a  compara¬ 
tive  theoretical  analysis,  including  prediction  of  amplitudes  using  a  very  simple, 
highly  idealized  linear  model.  The  experiments  were  performed  in  a  small 
laboratory  tank  in  which  a  densitygradient  was  created  by  maintaining  a  vortical 
temperature  gradient  in  water  and  the  mixed  region  was  formed  by  rapidly 
stirring  a  small  part  of  tho  fluid.  Tho  internal  waves  were  measured  by  a  ther¬ 
mistor  placed  near  ono  wall  of  the  tank.  Results  from  four  density  profiles  are 
illustrated.  Tho  theoretical  model  is  based  on  tho  theory  of  linear  modes  and  is 
used  to  predict  the  subsequent  motion  duo  to  tho  internal  waves  (amplitude  and 
phase)  in  terms  of  tho  maximum  turbulent  expansion  of  tho  mixed  region.  The 
forcing  function  in  this  case  is  composed  of  tho  initial  buoyancy  defect  moments 
in  tho  mixed  region  and  is  calculated  following  a  schemo  originally  proposed  by 
Lighthill  (1964,  private  communication).  Viscous  effects  are  included. 

Since  linearized  internal  wave  models  are  readily  amenable  to  mathematical 
treatment  it  is  believed  that  this  comparison  of  theory  and  experiment  is  doubly 
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valuable -it  notonly  provides&diackground  on  which  to  understand  the  experi¬ 
mental  results  but  also  . prorives  a  useful  indication  of  how  detailed  a  mathe¬ 
matical;  model  is  needed;  in  crderto  make  useful  quantitative  predictions  of 
resulting* wave  fields. 


2,  Discussion  of  erf  rimental  conditions  and  results 
Figure  1  (plate  l)  <G  ft.picture  of  the  'two-dimensional’  transparent  cell  that 
was  used,  -with  dim*.  .nsfons  30cm  wide,  7*3cm  deep  and  2-5 cm  thick  and  grid 
lines'2  cm  apart,  is  r  vhe  experiments  the  cell  was  completely  filled.with  distilled 
water  between  tliet  >tf  am  and  top  copper  strips.  Stable  stratification  of  the  water 
(more  dense  below  r  ,  m,  above)  was  produced  by  cooling  ihe  lower  copper  strip 
and  heating  the  r  one.  The  amount  of  cooling  and  heating  was  controlled  by 
the  polarity  and  , mount  of  electrical  direct  current  applied  to  commercial 
thermo-electric  -  /ts  (modern  Peltier  effect  devices)  attached  to  the  copper 
stiips.  For  effici-n.  andstable  operation  the  unattached  surfaces  of  the  M  cm 
thick.thermo-eb  vie  devices  were  held  at  a  constant  temperature  by  circulating 
water  through/,  ckots  as  shown  in-the  upper  part  of  figure  1. 


Temperatures  referred  to  4  cm  depth  (°C) 

Fiouue  2.  Toinpomturo  vs.  depth  profiles. - ,  profile  1; - ,  profile1 3 ; 

- ,  profilo  4;  -  ,  profilo  0. 

Figure  2  shows  four  different  profiles  of  temperature  vs.  depth  cstab'  .died  by 
using  various  combinations  of  electrical  current  through  tho  thermo-electric 
devices  (not  all  shown  in  figure  1).  For  example,  profilo  2  was  estaldrned  by 
passing  10  A  with  polarity  to  cool  tho  bottom  and  10  A  with  polarity  to  iieat  tho 
top.  Profile  3  required  5  A  cooling  current  on  the  bottom  and  5  A  heating  on  top. 
Profile  0  was  established  with  zero  bottom  current  and  2*5  A  heating  at'  the  top. 
Profilo  4  required  1 0  A  bottom  cooling  and  zero  top  current.  In  this  last  case  there 
was  slight  cooling  at  the  top  because  tho  constant-temperature  water  circulating 
through  tho  upper  water  jackets  was  somewhat  cooler  than  the  ambient,  tempera¬ 
ture.  The  cooling  was  by  conduction  through  the  structure  of  tho  thermo¬ 
electric  devices  attached  to  tho  upper  copper  strip.  It  took  about  2  h  after  tho 
start  of  current  flow  for  tho  temperature  profiles  to  stabilize  at  tho  values  shown 
in  figure  2.  (Profiles  2  and  5  are  omitted  from  this  report.) 
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Figure  3  shows  the  four  temperature  profiles  of  figure  2  converted  to  four 
corresponding  density  profiles.  The  conversion  was  made  using-handbook  data 
relating  pure  water  temperature  with  density,  together  with  a  reference  tempera¬ 
ture  for  each  temperature  profile  curve.  The  horizontal  scale  is  in  the  <rt  density 
unit  commonly  used  by  oceanographers  (Sverdrup  et  ai.  1942).  In  this  case 
crt  =  (p- 1)  x  103,  where  p  is  the  water  density  in  g/om3.  The  negative  abscissa 
values  are  density  expressed  in  parts  per  thousand,  and  density  increases  from 
left  to  right. 


Figure  3.  Density  (cr()  vs.  dopth  profilo3.  O — O,  profilo  1 ; 
x  -  -  x ,  profilo  3;  • — 9,  profilo  4;  ■  —|,  profilo  0. 

For  density  profiles  1  and  3  the  density  gradient  with  respect  to  water  depth, 
dcrtldz,  is  greater  at  the  top  and  decreases  with  depth.  For  profilo  6  the  density 
gradient  is  very  nearly  constant.  Profilo  4  shows  a  region  where  density  decreases 
with  dopth  near  the  top,  then  a  region  of  approximately  zero  gradient  which 
changes  to  aregion  of  approximately  constant-density  gradient  near  the  bottom. 

A  prominent  feature  near  the  centre  of  figure  1  is  a  1*3  cm  diameter  device 
which  is  a  non-rotatablo  mixer.  A  lover  system  connects  it  with  a  drive  mechanism 
that  will  move  the  mixer  forward  and  backward  about  ±1  cm  ten  times  in  about 
2-2-5  sec  on.  demand.  This  vigorous  back-and-forth  movement  will  generate  a 
pulse  of  turbulonco  which  will  first  expand  and  then  collapse  vertically  (Wu 
1969;  Schooloy  1908).  The  restoring  force  is  gravity  acting  on  p  and  dajdz .  It 
is  the  resulting  ohango  in  the  temperature  structure  with  time  after  mixing 
which  makes  it  comparatively  simple  to  study  the  internal  wave  structure  for  the 
various  temperature  or  density  profiles  of  figures  2  and  3, 

Figure  1  shows  a  small  thermistor  bead  projecting  about  2-5  mm  tlirougli  a 
small  hole  in  tho  upper  copper  strip,  directly  above  the  mixer.  The  thermistor 
was  connected  to  a  calibrated  recording  system  whioh  yielded  a  time  record  of  tho 
temperature  at  this  point  boforo,  during  and  after  mixing.  Since  tho  thermistor 
was  used  in  a  manner  which  was  relatively  insonsitive  to  velocity  it  recorded 
temperature  changes  duo  to  vertical  motion  of  tho  stratified  water  at  tho  point  of 
measurement.  A  second  thermistor  bead  is  also  shown  projecting  through  a  hole 
in  tho  copper  strip  near  tho  centre  of  tho  right  half  of  figure  1 .  Although  not  visible, 
there  are  small  holes  for  inserting  thermistors  at  four  different  places  along  tho 
copper  strip  on  eaoh  side  of  tho  centre-lino  (£). 
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Figure  4  shows  the  results  of  experiments.  In  figure  4  (a)  the  time  scale  proceeds 
•downward  for  40  sec  after  the  start  of  mixing  and  the  horizontal  scale  measures 
distance  from  zero  at  the  centre  of  the  upper  copper  strip  (<£)  arid  proceeds  16  cm 
to  the  fight.  The  curves  represent  smoothed  processed  data  from' several  experi¬ 
ments  using  profile  1.  The  long-dashed  curves  are  regions  where  the  thermistor 
recording  system  showed  positive  maximum  readings,  compared  to  the  situation 
before  mixing.  The  short-dashed  curves  show  regions  of  negative  maximum' 
temperature  readings.  The  solid  lines  aro  the  region*  , where  the  temperature  was 
the  samo  as  before  mixing. 


Distahco  right  of  centi-o-lino  (cm) 


Fiourk  4.  Experimentally  dotorminod  constant  phase  lines,  (a)  Temperature  profilo  1, 
(6)  tompprfttnfa  p.ofilo  3,  (c)  tomporaturo  profilo  0,  (d)  temperature  profilo  4. 

The  pattern  revealed  by  figure  4  (a)  is  that  of  a  series  of  internal  waves  initiated 
by  a  pulso  of  turbulence  caused  by  the  1*3  om  diameter  mixer,  4*6  cm  below  tho 
upper  water  boundary  and  on  tho  centre-line  of  tho  experimental  cell.  To  save 
space,  only  the  right  half  of  the  internal  wave  pattern  is  shown.  A  similar  sym¬ 
metrical  set  of  internal  waves  also  occurs  to  the  loft.  At  various  places  in  the 
negative  and  positive  regions  of  the  internal  waves  tho  approxi  mate  amplitudes 
of  tho  temperature  deviations  are  shown.  Since  tho  water  near  tho  top  was  warmer 
than  that  below,  a  negative  temperature  indicates  an  upward  component  of 
water  flow  and  a  positive  temperature  means  a  downward  component  of  flow. 
Tho  first  convergence  duo  to  tho  vortical  collapse  caused  tho  strong  positive  wave 
which  grow  from  about  +  0*  1 5  °C  noar  the  contro  to  -)•  0*5  °C  about  7*6  cm  from 
tho  centre. 

Internal  wave  absorbers  were  used  at  tho  right  and  loft  ends  of  tho  coll  (not 
shown  in  figure  1)  in  an  attempt  to  minimize  end  reflexions.  Howover  absorption 
was  not  complete  and  quantitative  data  was  not  attempted  when  it  was  thought 
that  reflexion  interference  might  bo  present. 

Figuro  4(6)  is  tho  samo  as  figure  4(a)  except  that  it  is  for  profile  3  instead  of 
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profile  1.  In  this  case  the  figure  shows  experimental  points  to  give  an  example  of 
the  amount  and  consistency  of  the  original  data.  The  temperature  gradient  is 
less  in  figure  4  ( 6 )  than  in  figure  4  (a)  diid  for  this  reason  the  initial  negative  wave 
was  too  small  to  be  measured.  The  first  wave  in  figure  4(6)  is  due  to  the  first  con¬ 
vergence  caused  by  the  vertical  collapse  of  the  initial  turbulent  divergence. 
The  internal  waves  for  the  profile  3  condition  (figure  46)  are  shown  to  travel  more 
slowly  than  for  the  stronger  temperature  gradient  profile  1  (figure  4o); 

Figure  4(c)*  is  for  the  weak  and  quite  linear  profile  6.  A  weak  temperature 
gradient  makes  temperature  variations  duo  to  the  internal  waves  more  difficult 
to  measure.  However,  enough  data  are  available  to  show  that  the  internal  wave 
speed  is  still  less  than  for  the  previously  discussed  stronger  gradients. 

Figure  4(d)  shows  the  internal  waves  when  profile  4  was  used.  In  this  case 
there  was  a  rather  weak  gradient  at  the  mixer  depth  with  an  isopycnal  inter¬ 
mediate  region  above.  At  the  top  the  water  was  slightly  colder  above  than  below, 
but  apparently  not  enough  to  induce  convection.  The  results  of  the  experiment 
show  slowly-travelling  internal  waves  where  the  collapse  of  the  initial  turbulent 
pulse  induces  first  a  negative  and  then  a  positivo  temperature  wave.  This  is  the 
reverse  of  figures  4  (6)  and  4  (c)  and  is  duo  to  the  change  in  the  sign  of  the  gradient 
near  the  top  for  profile  4  compared  to  the  case  for  profiles  3  and' 6. 


3.  Theoretical  analysis 

In  this  section  a  theoretical  model  which  is  based  on  the  conditions  which' 
prevailed  during  the  experimentation  is  discussed.  Calculations  using  this  model 
have  been  mado  and  are  presented  for  comparison  with  the  recorded  data. 

A  number  of  simplifying  and  idealizing  assumptions  have  been  used,  most  of 
which  are  common  to  the  first-order  study  of  internal  waves  (e.g.  Mowbray  & 
Rarity  1907).  Besides  the  usual  assumptions  of  linearity,  Boussinesq  approxima¬ 
tion  and  no  molecular  diffusion  of  those  properties  that  define  the  density,  the 
following  more  specialized  assumptions  aro  mado:  (i)  Molecular  viscosity  is 
included  by  means  of  a  perturbation  expansion,  following  Boro  (1968).  This 
assumption  will  bo  dealt  with  moro  fully  later,  (ii)  The  initial  disturbance,  and 
thus  the  resulting  flow  pattern,  is  strictly  two-dimensional,  (iii)  The  initial  velocity 
str ucturo  is  zero,  (iv)  The  end  walls  aro  infinitely  distant,  (v)  Tlw analysis,  later  on, 
is  particularized  to  the  case  in  which  the  Brunt-Vaisala  frequency  is  independent 
of  depth.  It  is  unfortunate  that  more  information  was  not  available  about  tho 
degreo  of  mixing  prior  to  collapse.  Howovcr,  it  is  oxpcctcd  that,  at  least  over  the 
area  of  tho  stirrer,  mixing  was  virtually  complete  and  any  incompleteness  was 
near  tho  edgo  of  tho  mixed  region.  If  this  wore  tho  caso,  it «  juld  show  up  in  tho 
resulting  wave  field  mainly  as  a  deficiency  of  high-order  modes.  It  will  be  shown 
below  that  lor  tho  comparison  with  tho  experimental  data,  low-o^dor  modes  aro 
dominant  (owing  to  geometrical  coupling  and  viscous  effects).  For  the  present 
problem  therefore,  the  simplest  mathematical  representation  of  tho  initial  state 
will  be  used,  namely  a  perfectly  mixed  circular  region , 

Of  all  these  assumptions,  (iii)  and  (v)  are  probably  vho  most  unlike  tho  experi¬ 
mental  conditions  for  short-term  observations.  Initially,  and  very  near  tho 
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*  edge  ’  of  the  mixed  region,  all  assumptions  (except. possibly  (iv))  are  expected  to 
break  down.  However,  this  breakdown  should  apply  mainly  to  very  short  length- 
scales,  and  since  their  group  velocities  ate  low  and  viscous  damping  high  they  are 
not  expected  to  contribute  materially  to  the  flow  pattern.  Also,  the  turbulent 
velocity  structure  existing  initially  will  be  most  energetic  in  length-scales  less 
than, the  initial  diameter  of  the  mixed  region  and  will  be  almost  completely 
contained  within  the  mixed  region;  thus  the  direct  coupling  between  this  struc¬ 
ture  and  the  resulting  wave  field  vail  be  low  (especially  for  waves  long  compared 
with  the  initial  diameter). 

Non-linearities  in  the  resulting  wave  field  may  be  important  for  a  short  time 
after  the  beginning  of  collapse  and  near  the  mixed  region.  An  upper  estimato  for 
the  vertical  amplitude  of  wave  motion  is  (Wu  1969),  where  a  is  the  radius  of 
the  mixed!  'gion,  thus  the  ratio  of  the  non-linear  terms  in  the  equations  of  motion 
to  a  typical  linear  term  is  given  by  \ak,  where  k  is  the  wave-number  of  the  wave- 
field.  Prom  figure  4  and  table  1  the  largest  value  of  this  parameter  is  0-35  (profile  1). 
Pur  regions  or  times  removed  from  the  onset  of  collapso  this  (parameter  is 
expected  to  be  substantially  reduced  (figure  7  (a)  or  Wu  1969). 

The  position  at  whicVthe  measurements  were  taken  is  very  near  a  solid  boun¬ 
dary  so  there  is  the  possibility  of  a  substantial  modification  of  the  inviscid  wave- 
pattern  by  boundary-layer  effects'.  Also,  since  the  Brunt-Vaisala  frequency  is  of 
order  1  per  sec  there  is  tho  likelihood  of  apprecitii  le  viscous  dissipation  within  the 
body  of  the  fluid  and  by  the  walls.  It  was  not  expedient  to,  perform  a  new  set  of 
experiments  (the  ones  reported  hero  were  done  some  years  before  tho  theory  was 
undertaken),  so  instead,  these  viscous  effects  are  included  in  tho  theoretical  model. 
To  accomplish  this,  the  perturbation  scheme  given  by  Doro  will  bo  used.  It  is 
directly  applicable  to  the  present  problem  and  so  only  the  results  of  that  investi¬ 
gation  will  bo  used.  It  is  thus  necessary  to  solve  only  the  inviscid  equations  and 
apply  tho  viscous  corrections  lato;-. 

With  the  previously  mentioned  assumptions,  tho  inviscid  equations  of  motion 


du  8w 

(1) 

to+¥  =  °’ 

.dll  dp 

',»(z)8i+5J‘0' 

(2) 

(3) 

(•1) 

with  x  horizontal  and  perpendicular  to  tho  axis  of  tho  mixed  region,  s  positive 
downwards, p0(z)  the  density  profilo  before  mixing,  p  tho  instantaneous  difference 
between  the  total  density  and  p0,  p  tho  instantaneous  difference  between  t..e 
total  pressure  and  the  hydrostatic  head,  g  tho  acceleration  of  gravity,  u,  w  the 
velocity  components  in  tho  x,  z  plane,  and  t  tho  time.  Tho  boundary  conditions  aro 

w  =  0  at  z  =  0,D  for  all  time,  (5) 

w  remains  bounded  as  aJ-»±oo.  (6) 
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The  initial  conditions  are 

p,u,w  =  0  for  all  x,z  (7) 

and  p  =  p0(z)f(x,z).  (8) 


The  problem  is  amenable  to  treatment  by  a  one-sided  Fourier  transform  in 
time.  Thus,  if 

w=j  e~iuw(t)dt  (9) 

(and  similar  representations  for  all  the  other  dependent'  variates),  equations  (1) 
to  (4)  and  (8)  reduce  to  the  following  equation  (using  the  B outlines q  approxi¬ 
mation):  a*®  /  me*w_  gd*f 

dz2  +  \  a 'j*Jdx2~  «*&*’  {) 


where  S(H  -  (£§ f 

is  the  Brunt-Vaisala  frequency. 

To  reduce  the  problem  further,  lot 

S^„(2)Fn(*), 

n 


(ii) 


where  $Jn(3)is  an  eigenfunction  from  an  orthonormal  set  defined  by 

Wi+a"(S“1)^n=:()  ^n(z)  *  0  at  s  =  0,1))  (12) 

d2F  a  CD  d2f  d?I 

,„dthm  (13) 


Since  1  is  a  direct  representation  of  the  buoyancy  at  time  zero  it  must  vanish 
outside  the  initial  mixed  region.  It  is  also  necessary  that  <j>n(z)  be  the  nth  member 
of  a  complete  sot;  otherwise  an  integral  term  (representing  a  continuous  dis¬ 
tribution  of  modes)  is  required  on  the  right-hand  side  of  (11).  It  is  shown  by 
Couraut  &  Hilbert  (1963)  that  <j>n(z)  is  complete  if  w  possesses  continuous  first 
and  piecewise  continuous  second  derivations  (in  z),  i.o.  from  (10),  if  N  (or/)  is  at 
worst  piecewise  continuous  in  z.  Thus  the  solution  of  (13)  satisfying  (6)  is 

Fn(x)  —  I  -  \ian  ^eianrJ  (y)  dy + e~ia^J  eianvj(?/)  d?/J .  (14) 

Let  the  initial  mixed  region  bo  symmetric  in  x  with  its  midpoint  at  (0, 30)  and  its 
maximum  horizontal  extent  given  by  2a.  Then  the  total  wave  field  exterior  to  the 
initial  mixed  region  can  bo  obtained  from 

Fn{x  •?  f  cos  (ani/)I(ij)dt]  (|a:|>a).  (16) 

Using  the  invorso  transform  to  (9), 

W  =  5rf„  (“,‘anc“fa"w+<wt5in(s:)  [o  003 do),  (10) 

where  e  is  chosen  to  ensure  that  the  integration  is  below  the  singularities  of 
the  integrand. 
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Internal  loaves  in  a  im-dimensional  mixed-region 
The  initial  buoyancy  function^)  is  given  by 

I  =  ~ij^f{x,z)<f>n{z)dz, 

where,  for  a  circular,  completely  mixed  region, 

j{x,  z)  -  (PoizJlPoiz)  - 1 )  H(x + a)  H(a  -  x)  U(z  -  z0 + («2  ~  #  ( V* + («2  -  ^)h)._ 

(17) 


Figure  6.  Thoorotioal  modol  of  tho  initial  density  vs.  dopth 
through  tho  oentro-lino  of  tho  mixod  rogion. 


Hero  Zq  is  the  depth  of  tho  centre  of  the  region,  H  is  the  Heaviside  unit  function 
and  po(zi)  is  tho  density  of  tho  mixed  region  (see  figure  6). 

In  terms  of  tho  Brunt-Vfiisula  frequency1, 

or,  using  the  Boussincsq  approximation, 


Thus, 


1  = 


Po(z)  J 9 
ff(W)fl(q-3;) 


W{z')dz'dz. 


(18) 

(10) 


The  general  solution  represented  by  equations  (19),  (16)  and  (12)  will  now  bo 
restricted  to  tho  present  experimental  ease.  Tho  density  profiles  shown  in  figuro 
3  indicate  that  tho  approximation  IV2  =  constant  should  provide  most  of  tho 
pertinent  features  of  tho  internal  wavo  field.  With  N  =  N0  and  zx  =  z0, 


,  .  .  2 lain (nnzID) 

PnW  “  l))i- 

(20) 

nn 

“  D{NH(oa--l}i' 

(21) 
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J%os(«„,)7 'Mr  — f *)]' “(tt)  -4 (sjr^W) '  (22) 


where  J2  is  the  Bessel  function  of  second  order.  Therefore,  interchanging  the 
order  of  integration  and  summation  in  (18)  (which  is  permissible  by  the  rules 
governing  Fourier  series  representing  generalized  functions), 


w  =  xw 


„  ~  Bm(nj)zlD)cos[nnz0ID)  C+aj-<c  T 

“  >. - : - I  <X„t/2 

nn  J-oo-/c 


S 

n  -»1 


/  an7r  \ 


gfut-flzl  *ndo>. 

(23) 


Tliis  expression  can  be  simplified  by  contour  integration.  For  l  >  0  it  can  be 
jhown  that  there  is  no  modification  to  the  value  of  this  integral  if  the  integration 
path  is  closed  by  anrinfinito  somicirclo  in  the  upper  half  plane.  For  t  <  0  the 
lower  half  plane  can  bo  used.  Thus  the  value  of  w  arises  only  from  integration 
around’ the  singularities  io  —  ±N0  (in  an).  Since  these  are  branch  points  a  cut  is 
necessary  and  it  is  convenient  to  take  the  cut  along  tho  real  o  axis  between  +  AT0. 
Integration  around  tho  branchpoints  produces  no  contribution,  so  (23)  reduces  to 
an  integration  along  the  bottom  of  tho  cut  from  -N0  to  +N0  and  an  integration 
along  the  top  in  the  opposite  direction.  In  appendix  A  it  is  shown  that 


d{a\)ld{wr)  >  0 

near  the  real  o>  axis  and  that  this  condition  together  with  the  condition  of 
boundedness  at  large  j«J  leads  to  the  restriction  t'.iat  the  real  parts  of  w  and  an 
must  havo  the  same  sign  for  (o  slightly  below  the  real  axis.  Therefore,  on  the 
bottom  of  the  out  <  0  for  (v%  0  and  on  top  an  $  0  for  (o%  0. 

Finally,  with  e  »  0, 


xv  -  — 


•2T^,3,n( 


117TZ\ 

W  003 


ms, 


X'J4awrlD{\-(o-lNl)h) 


x  sin  (ol  cos 


(Ag/&>2-  1)1 
/  nnx  \ 
\D(NUcv> -l}i) 


do.  (24) 


For  l  <  0  there  are  no  singularities.  Thus 


tiis  0  for  1  <  0. 


(25) 


Using  tho  linearity  assumption  and  tho  transformation  w  =  AT0cosd,f  the 
amplitudo  f  (= jxvdt )  is  given  by 


in 


r!*. 

/  amt  \ 

Ur  l  v)j 

I.  • 

'•\D  sin  O) 

(cos 


(  mix  \ 

m^ndj 


Tliis  represents  t  he  inviscid  solution. 


dO, 

(20) 


t  It  van  cosily  bo  seen  that  to  tins  approximation  0  represents  the  angle  between  the 
vortical  and  tho  local  lines  of  constant  phase  or  ‘crest*  linos  of  that  particular  frequency 
component.  If  0  ss  0  the  lines  of  constant  phase  are  vertical  and  if  0  =  \tt  the  lines  are 
horizontal. 


Wav©  height  (10~3  cm) 
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Internal  waves  in  a  two-dimensional  mixed  region 

To  incorporate  the  most  important  effects  of  viscosity  it'is  necessary  to  return 
to  equation  (16)  and  perturb  the  frequency  o)  in  eiut  and  to  modify  £  to  allow  for 
boundary-layer  effects.  The  resulting  equations  are  obtained  in  appendix  B. 
Some  general  characteristics  of  tlio  form  that  tho  internal  wave  pattern  possesses 
will  now  be  outlined. 

The  finite  depth  of  the  model  1ms  resolved  the  pattern  into  an  infinite  series  of 
modes.  (This  situation  remains  true  for  an  infinite  depth  if  N  is  piecewise  con¬ 
tinuous  and->  0  as  2->co).  The  horizontal  inviscid  group  velocity  of  the  nth  mode, 
Cgn,  is  given  by 


0  5  l.n  15  0  2  4  6  8 


(a)  Distanco  in  units  of  D  (6)  Timo  in  units  of  tho  Brunt-Vuisiilti 

period 

Figuhe  0.  Illustration  of  tho  mode  structure  of  an  inviscid  internal  wave  Hold.  Tho  inodes 
lmvo  been  separated  to  show  comparative  fa  lures,  (a)  wuvo  h^Sit  at  an  instant  as  a 
function  of  distance  from  tho  centre-line,  N^t  =•  50j  (6)  wavo  height  at  ono  particular  point 
as  a  function  of  t  imo  from  beginning  of  tho  collapse,  x/D  =  2.  (For  thisoxnmplo;  AT0  a  0*25 
seo*1,  a  =  2-05  cm,  D  =  7-3  cm,  z  =  0-3  om,  z0  ss  4*5  cm  and  viscosity  =  0.) 

Also,  the  horizontal  phase  velocity  0H  is 

~  (0  NnD  .  n 

(?  =  —  =  — 2—sm  0.  (276) 

n  an  717T  v  ' 

Tho  maximum  group  velocity  N0Djn  (which  is  tho  same  as  the  maximum  phase 
velocity)  increases  directly  with  A70.  This  is  in  accordance  with  the  experimental 
observation  that  tho  wavo  velocity  increases  as  the  density  gradient  (or  jV0) 
increases.  Also,  the  horizontal  extent  of  the  wave-pattern  fur  each  mode  is  limited 
at  any  instant  and  at  higher  mode  numbers  the  limitation  appears  nearer  the 
centre-line  of  the  mixed  region.  Within  this  horizontal  limitation  a  wave-field 
stretching  back  to  the  centre-line  can  exist. 

A.  stationary  phase  reduction  of  equation  (26)  (with  N0  large)  indicates  that  a 


r:  i  r  ' "  v  - ; 
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stationary  phase  contribution  exists  for  every  point  within  this  possible  wave- 
field.  Equation  (27)  gives  the  value  of  0  at  which  the  phase  of  equation  (20) 
is  stationary  for  \x\\l  =  Cgn.  As  (z|  ->  0,  0->  0,  the  horizontal  wave-number 
M7r/Z>  tan  oo  and  to->N0.  Therefore  each  mode  gives  rise  not  to  a  localized 
group  of  waves  but  to  a  field  extending  from  the  centre-line  to  the  limiting 
extent  for  that  particular  time.  The  wavelength  within  the  field  increases  from 
zero  to  co  over  this  same  region  Eigure  0  (a)  illustrates  this  behaviour. 

As  a  function  oftime  at  a  given  position  the  wave-pattern  is  essentially  zero 
until  the  most  rapid  (n  =  1)  mode  arrives  (at  t,  say)  and  exhibits  a  continuing 
oscillation  in  time  with  a  frequency  rising  from  zero  asymptotically  at  N0.  At  2 1 
the  second  mode  arrives  and  exhibits  a  similar  behaviour  in  frequency  and  so  on 
for  all  the  modes.  This  is  shown  in  figure  6(b). 

If  only  one  mode  is  present  the  wave-field  will  exhibit  dispersion  characteristics 
pertaining  to  that  mode.  If  more  than  one  mode  is  dominant,  as  is  true  for  the 
example  used  in  figure  0,  interference  cfTccts  will  also  be  apparent. 

The  damping  effects  of  viscosity  are  most  pronounced  for  high  modes  and  high 
frequencies  (0  «  0).  Waves  nearest  the  centre-lino  at  any  given  instant  or  waves 
appearing  at  later  times  at  any  given  position  are  reduced  the  most.  The  local 
frequency  is  also  reduced  by  viscosity. 

All  these  general  characteristics  aro  in  agreement  with  the  measured  patterns 
shown  in  figure  4  oven  though  these  do  not  satisfy'  the  condition  Ar0t-> co.  The 
only  major  excopffions  occur  in  tho  vicinity'  of  the  ini  tial  mixed  region.  In  figure 
4(a)  the  distance  between  zeros  along  tho  10  sec  line  increases  away'  from  the 
centro-linc.  Tho  same  behaviour  can  bo  seen  in  4  (b)  along  tho  20  sec  line.  Also,  the 
time  between  zeros  along  the  5  cm  line  in  figure  4(b)  decreases  towards  later 
times  (although  in  figure  4(a)  the  bottom  curves  indicate  a  slight  opposite 
tendency).  However,  not  only’  may  asymptotio  arguments  be  expected  to  fail 
there  and  viscous  offccts  dominate,  but  the  region  |$|  <  a  (for  all  z)  has  been 
previously’  excluded  from  the  domain  of  applicability'  of  the  solution.  It  is  also 
apparent  from  tho  measurements  that  a  simple  wave  field  exists,  implying  the 
presence  of  only  a  few  m*  des. 

To  accomplish  a  direct  comparison  of  tho  theory  with  experimental  results 
the  temperature  excursion  from  ambient  was  calculated  using  AT  =  £(< iTjds ) 
evaluated  at  z,  the  depth  of  the  thermistor  bead,  for  each  profilo  shown.  Tho 
values  used  in  equation  (B3)  are  listed  in  table  l.  Each  of  these  values  was  ob¬ 
tained  from  the  experimental  data.  The  listed  values  of  A'0  are  the  averages  of  t  he 
Brunt- Vaisfila  frequency  over  the  total  depth  of  the  tank  for  each  profile.  In 
eaeh  ease  the  average  was  calculated  from  an  exponential  curve  of  N(z)  fitted 
with  a  least-square  technique  to  individual  values  of  N  obtained  from  figure  3. 
The  accuracy  of  N0  is  estimated  at  +  lo%.  Maximum  vertical  expansion  of  the 
turbulent  pulses  was  determined  by  dye  and  cinematography  techniques 
(Schooley’  1008).  One-half  the  expansion  determined  a  to  an  accuracy  of  about 
±10%. 

The  values  of  viscosity’  listed  were  obtained  from  standard  handbook  tables 
and  pertain  to  the  measured  values  of  temperatures:  i*  corresponds  to  an  average 
throughout  the  tank,  r  corresponds  to  the  value  at  the  thermistor  depth. 
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It  was  not  expected  that  mode  1  would  be  important  since  it  represents  fluid 
travelling  all  up  or  all  down  in  a  given  column.  Instead  the  initial  buoyancy 
pattern  is  such  that  it  would  couple  well  into  modes  which  represent  a  downward 
velocity  above  z0  and  an  upward  velocity  below  z0.  If  z  were  exactly  half  of  D,  the 
flow  would  therefore  be  mode  2  with  contributions  from  the  other  oven  ordered 
modes  because  of  the  initially  round  shape.  In  fact,  in  (2G)  all  odd  modes  do 


Profilo 

a  (cm) 

N0  (sco-1) 

dT/dz 

(°C/cm) 

l'(cgs) 

z  (cm) 

1 

2*00 

HO 

11*1 

f*000 

0*250 

3 

2*05 

0*08 

3*8). 

0*001) 

0*300 

0 

2*35 

0*41 

1*18 

0*010 

0*325 

Tabus  1. 

-o  =  4*5  cm,  D 

=  7*3  cm, 

v  ss  0*01  cin/soo4 

vanish  under  this  condition.  With  tho  present  ratio  of  zJD  =  0*610  it  is  expected 
that  mode  3  will  also  bo  important  since  the  cigonfunction  for  n  **  3  has  a  zero 
at  zjD  -■  §.  A  large  contribution  from  modes  2  and  3  is  apparent  in  the  example 
shown  in  figure  0.  For  tho  calculations  used  in  tho  comparison  with  tho  m  aaured 
data  enough  modes  were  included  in  the  summation  to  ensure  stability  of  one  part 
in  103  in  the  final  valuo  of  £.  The  necessary  number  of  modes  ranged  from  4  at. 
largo  |x|  and  £  to  30  at  small  |a:|  and  £.  Typically,  mode  2  was  dominant  with  some 
substantial  contributions  from  modes  3  and  5. 
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Figure  7.  Temperature  excursion  vs.  time.  theoretical  curves; experimental  vnluiv. 
Tho  position  of  eaoh  lino  is  in  scale  with  its  distance  from  tho  centre-lino.  The  distances  in 
centimotres  are  shown  abovo  each  lino. 


Tho  comparison  is  shown  in  figure  7.  The  dotted  lines  represent  measured 
data,  the  solid  lines  are  the  theoretical  curves.  At  least  two  sots  of  measured  data 
were  obtained  for  each  profile.  Tho  time  origin  for  the  calculated  curves  has  had 
one  adjustment  for  each  profile  to  provido  a  visual  ‘best’  fit  to  the  measured 
curves.  This  is  considered  permissible  becauso  the  instant  when  collapse  began  is 
not  known  experimentally.  This  is  becauso  the  concept  itself  is  an  idealiza¬ 
tion:  collapse  can  occur  during  tho  entire  stirring  and  turbulent  expansion 
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interval.  Also,  a  small  change  in  N0  or  ly  variation,  in  the  Brunt-Vaisula 
froquonoy  with  dopth  will  affeot  tho  velocity  of  the  waves -without  appreciably 
affooting  their  amplitudes,  thus  effectively  producing  a  difference  between  experi¬ 
mental  and  theoretical  time  origins.  Tho  time  shifts  actually  used  are  1*0, 1*0  and 
4*0  see  for  profiles  1, 3  and  0  respectively.  As  in  figure  4,  the  tin  o  origin  refers  to 
tho  instant  tho  mixer  was  turned  on.  For  each  curvo  tho  distanco  in  centimetres 
from  tho  contro-lino  to  tho  buso-lino  at  which  thoso  measurements  were  obtained 
is  shown  by  tiro  number  at  tho  top  of  each  figure.  Tho  excursions  from  the  baso- 
lino  ropresonttho  temperature  deviations  from  ambient  and  aro  scaled  at  the 
bottom  of  each  figure.  For  profilo  1  tho  fit  is  good  only  for  tho  first  wavo.  After 
that  the  theoretical  curves  predict  moro  waves  than  aro  experimentally  found 


Fniuw:  S.  Calculated  tempornluro  excursions  for  profilo  3  with  zero  viscosity.  Tlio  number 
above  each  lino  is  its  horizontal  distanco  from  tho  contro-lino  in  centimetres. 


(oxcopt  for  tho  measurements  at  12-7  cm).  Nevertheless  tho  amplitudes  of  tho 
theoretical  excursions  aro  very  similar  to  the  measurements.  For  profiles  3  and  0 
tho  measurements  and  tho  theoretical  ettrves  aro  in  reasonably  good  agreement. 

Tho  invisoid  solution  for  profilo  3  is  shown  in  figure  8,  It  can  bo  scon  that  for 
this  scalo  of  motion  viscosity  is  indeed  important  and  over  tho  recorded  time 
interval  reduces  tho  amplitude  of  tho  wave  typically  by  a  factor  of  two.  Tho 
degree  of  matching  between  (B  1)  and  (22)  was  tested  by  sample  calculations  on 
profiles  1  and  0.  Equation  (B  1)  was  used  in  (B  3)  with  i»  =  0  (i>  0)  nnd  compared 

with  tho  results  for  $n  in  place  of  in  (B  3).  Both  functions  gavo  the  same  results 
with  differences  of  about  5  %. 

Equation  (20)  indicates  that  tho  wavo  amplitude  is  fairly  strongly  dependont 
on  tho  radius  of  tho  initial  mixed  region.  If  «//)  1, 

JtiannjDmiO)  ss  (annlDthinO)*  x  £ 
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for  low  mode  numbers  and  low  frequencies.  In  this  case  foe  a*.  If  the  mode  number 
of  the  dominant  waves  is  largo  enough  or  the  frequency  of  the  dominant  waves  is 
small  enough,  or  if  ajD  =  0(1),  J^ann/D  sin  0)  is  essentially  independent  of  a  and 
so  foe  a2. 


4.  Conclusions 

It  has  been  shown  that  the  internal  wave  amplitudes  generated  by  the  two- 
dimonsional  collapse  of  an  initially  turbulent  rcg'on  can  bo  predicted  (at  least 
near  a  'boundary)  using  a  highly  idealized  linear  model.  The  accuracy  of  the 
prediction  is  only  loosely  determined  but  for  those  measurements  it  is  approxi¬ 
mately  the  same  as  the  experimental  scatter.  The  discrepancy  in  time  of  arrival 
of  the  waves  botwcon  theory  and  oxporimont,  as  is  indicated  by  the  adjustment  in 
the  timo  origin  of  tho  theoretical  values,  is  beliovedto  bo  due  mainly  to  two  effects: 
variations  in  the  Brunt-Vfiisala  frequency  from  tho  assumed  constant  profile  and 
non-zero  mixing  timo. 


Appendix  A 

The  eigenvalue  an  and  the  eigenfunction  <fin(z)  are  defined  by 

$  +  a?‘(5‘)~1)^  =  0  " 0  at  Z  =  0>D)  (A1) 


and 


jre-)w*- 


i. 


(A  2) 


(Equation  (A  1)  is  tho  same  as  (12)  in  tho  main  part  of  the  paper.)  Differentiate 
(A  1)  by  (o-  and  set  8<^Jd<o-  =  i/rn, 

d2i/rn  .  (N-  \  .  dot-.  /A72  A,  o&N2  ,  n 

dz-  +a'V-  1)^"+dw‘-  (<e2  7  (o*  ^n_0,  (A3) 

Multiply  (A3)  by/0*.  the  complex  conjugate  of  0,„  and  integrate  over  all  z,  then 

Take  tho  complex  conjugate  of  (A  1),  multiply  by  ijrn  and  integrate  over  all  w. 
Since  (0  is  real  an  is  real, 


<A5) 


and,  from  (A  4), 


Therefore  (da2) Ido-  >  0  and,  in  fact,  since 
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(by  A  3),  which  is  greater  than  1,  da\!dw 2  >  a?,/<o2.  This  shows  that  the  horizontal 
group  velocity  is  always  less  than  the  horizontal  phase  velocity. 

Let  oj  =  (or-ie  and  an  =  anr+i£.  Then  for  small  e, 


However,  in  order  that  F(x)  remains  bounded  as  «->  +  oo,  f)  must  be  negative 
(equation\14).  Therefore,  since  e  >  0  (equation  16)  and  da\ldo)”  >  0,  <ur  and  «nr 
must  have  the  same  sign. 


Appendix  B‘. 

Using  the  perturbation  expansion  set  up  by  Doro,  equation  (20)  .represents  the 
inviscid  outer  solution  to  the  problem.  The  lowest  order  ‘inner  ’  solution,  near  the 
top  boundary,  is  given  by 


,  2l  fnnz 


+  ^(^)  *  p  ~  2f  cos  (z(w/2i>)i  +  £77) 

+  t[l  -  2i  cos  (z((oj2i>)i  -  $jr)]l ) , 


(Bl) 


where  v  is  the  local  kinematic  viscosity.  It  is  also  shown  by  Dorc  that  the  fre¬ 
quency  in  equation  (16)  must  be  replaced  by  a  perturbed  frequency  a)p.  To  second 
order  in  the  expansion 


(o„ 


(B  2) 


where  V  is  the  average  of  the  kinematic  viscosities  at  z  -  0  and  z  =  D. 

If  (B 1)  and  (B  2)  are  substituted  in  (23)  and  the  previously  outlined  steps  are 
taken,  an  equation  for  £  analogous  to  (26)  is  obtained.  In  a  e.mviified  form 


0  co 

£  =  =  S  cos 

Tf  n»l 


'nnzA  rX'On  .  / _ nnn  \ 

,d)} o 


c-«(i 


x  [^rcos  (o)rt)  cos  (anx) + jj^sin  (o>rt)  cos  (ctnx)\  dw,  (B  3) 
=  ^r-r  tfti  and  (i>p  =  o)r((o)  -  i<at(a>). 


where 


IMioio^mph  of  tin*  IjiiiU. 


o 
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